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We report on a search for gravitational wave bursts using data from the first science run of the LIGO detectors. 
Our search focuses on bursts with durations ranging from 4 ms to 100 ms, and with significant power in the LIGO 
sensitivity band of 150 to 3000 Hz. We bound the rate for such detected bursts at less than 1.6 events per day 
at 90% confidence level. This result is interpreted in terms of the detection efficiency for ad hoc waveforms 
(Gaussians and sine-Gaussians) as a function of their root-sum-square strain h rss ; typical sensitivities lie in the 
range hrss ~ 10~ 19 — 10~ 17 strain/\/Hz, depending on waveform. We discuss improvements in the search 
method that will be applied to future science data from LIGO and other gravitational wave detectors. 

PACS numbers: 04.80.Nn, 07.05.Kf, 95.30.Sf, 95.85.Sz 
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I. INTRODUCTION 

Gravitational wave bursts are expected to be produced from 
astrophysical sources such as stellar collapses, the inspirals 
and mergers of compact binary star systems, the generators 
of gamma ray bursts, and other energetic phenomena. Upper 
limits from searches for gravitational wave bursts with reso- 
nant bar detectors have recently been reported in [1—4], and 
results using interferometric detectors are published in [5, 6]. 
A new generation of detectors based on laser interferometry 
has been constructed, aiming for direct detection with broad- 
band sensitivity. These include the three LIGO detectors [7] 
described briefly in section II, as well as the British-German 
GEO 600 detector [8, 9], the Japanese TAMA 300 detector 
[10], and the French-Italian VIRGO detector [11], forming a 
worldwide network. In the summer of 2002, all three LIGO 
detectors were brought to their design optical configuration. 
After a series of engineering runs, the LIGO, GEO 600 and 
TAMA 300 detectors operated in coincident observation mode 
for the first time (Science Run 1, or SI) for two weeks in 
August-September 2002. 

Although the LIGO detectors were far from their design 
sensitivity, the quality of the data was sufficiently high to ex- 
ercise the first generation of analysis procedures for various 
types of gravitational wave searches, including searches for 
chirp gravitational waves from compact neutron-star binary 
inspirals [12], quasi-monochromatic gravitational waves from 
pulsar J 1939+2 134 [13], and broad-band stochastic back- 
ground gravitational radiation [14]. In all these analyses, a 
well-defined astrophysical model was assumed. In this paper 
we report on a search (using the LIGO detectors only) for un- 
modeled gravitational wave bursts that might originate from 
supernovae in our galaxy, mergers of binary stellar-mass sys- 
tems, gamma ray burst engines, or other energetic sources. 
The waveforms of gravitational waves from such sources are 
poorly known, so we employ data analysis algorithms which 
can, in principle, identify bursts with a broad range of possible 
waveforms. 

The first detection of gravitational wave bursts requires sta- 
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ble, well understood detectors; well-tested and robust data 
processing procedures; and clearly defined criteria for estab- 
lishing confidence that no signal is of terrestrial origin. None 
of these elements were firmly in place as we began this first 
LIGO science run; rather, this run provided the opportunity 
for us to understand our detectors better, exercise and hone our 
data processing procedures, and build confidence in our abil- 
ity to establish detection of gravitational wave bursts in future 
science runs. Therefore, the goal for this analysis is to produce 
an upper limit on the rate for gravitational wave bursts, even if 
a purely statistical procedure suggests the presence of a signal 
above background. It should also be noted that the sensitivi- 
ties of the three LIGO detectors during S 1 were several orders 
of magnitude worse than required for plausible detection of 
bursts from astrophysical sources such as supernovae in our 
Milky Way galaxy [15]. 

In this search we focus on short (4 ms to 100 ms) bursts 
in the LIGO sensitivity band (roughly 150 to 3000 Hz), with 
sufficiently high strain amplitude to be observed over the de- 
tector noise. We make no other assumptions about the nature 
or origin of the burst. We apply software algorithms to the 
LIGO detector data stream to detect such bursts. In order to 
suppress false signals from fluctuations of the detector noise 
we require temporal coincidence of detected burst events in 
all three LIGO detectors. We estimate the rate of accidental 
coincidences by studying the number of time-shifted coinci- 
dent burst events, and look for a statistically significant excess 
of coincident burst events at zero time shift. In light of the 
discussion in the previous paragraph, our goal for the search 
presented here is to set an upper limit on the rate of excess co- 
incident bursts, given the detectors' level of sensitivity during 
the S 1 run. 

In order to interpret our upper limit on the rate of burst 
events, we evaluate the efficiency of our search algorithms 
for the detection of simulated bursts injected into the data 
streams, using simple, well-defined waveforms (Gaussians 
and sine-Gaussians). We obtain curves of triple-coincidence 
detection efficiency as a function of gravitational waveform 
peak amplitude at the Earth, averaged over source direc- 
tion and incident wave (linear) polarization. We then com- 
bine our gravitational wave burst rate limits with these effi- 
ciency curves, yielding rate-versus-strength regions that (for 
the waveforms that we have examined) are excluded at the 
90% confidence level or higher. 

The paper is organized as follows. In section II we briefly 
describe the LIGO detector array and the data obtained from 
the first science run, with emphasis on those characteristics 
most relevant for a search for short gravitational wave bursts. 
In section II B we briefly describe the S 1 run. In section II C 
we describe the data quality requirements that were applied 
to the S 1 data sample, and present the subset of the data used 
for this search. In section III we describe our data processing 
pipeline, including the event trigger generation, event vetoes, 
and the time coincidence requirement. We present the results 
of two independent pipelines, based on the burst detection al- 
gorithms discussed in section IIIC. In section IV we estimate 
the background (accidental coincidence) event rate. In section 
V we evaluate the efficiency for the detection of bursts mod- 
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eled with simple ad hoc waveforms, and compare that with 
expectations. In section VI we present our limit on the ob- 
served excess event rate. We combine this with our efficiency 
curves as a function of signal strength, excluding regions in 
the rate versus signal strength plane. We also discuss the most 
significant systematic errors in these measurements. We sum- 
marize these results in section VII. Finally, we outline our 
plans to improve and expand our search methodology using 
data from subsequent observation runs. 



H. DETECTORS AND DATA SET 

A. The LIGO detectors 

All three LIGO detectors are orthogonal arm Michelson 
laser interferometers. The LIGO Hanford Observatory oper- 
ates two identically oriented interferometric detectors which 
share a common vacuum envelope: one having 4 km long 
measurement arms (referred to as HI), and one having 2 km 
long arms (H2). The LIGO Livingston Observatory operates a 
single 4 km long detector (LI). The two observatories are ap- 
proximately 3000 km apart, corresponding to 10 ms of light 
travel time. The detectors are approximately co-aligned, so 
that a gravitational wave should appear with comparable sig- 
nals at both sites. The principles underlying these laser inter- 
ferometer gravitational wave detectors are discussed in [16]. 
A more detailed description of the LIGO detectors can be 
found in [17]. 

These detectors aim to detect gravitational waves by inter- 
ferometrically monitoring the relative separation of mirrors 
which play the role of test masses, responding to space-time 
distortions induced by the waves as they traverse the detec- 
tors. The effect of a quadrupolar gravitational wave is to pro- 
duce a strain in space, impinging upon the detector and thus 
displacing the mirrors at the ends of the arms by an amount 
proportional to the arm length. For gravitational waves in- 
cident from directly overhead or below, and polarized along 
the arms of the detector, the mirrors at the ends of the two 
arms experience purely differential motion. Waves incident 
from non-optimal directions and/or polarizations can also in- 
duce differential motion; the "antenna pattern" is discussed in 
section V C. 

Each interferometer is illuminated with light from a 
Nd:YAG laser, operating at 1064 nm [18]. Before the light 
is launched into the interferometer, its frequency, amplitude 
and direction are all stabilized, using a combination of active 
and passive stabilization techniques [18, 19]. The light is sent 
through a beam splitter towards both arms. In each arm, a 
pair of mirrors (the "input test mass" and "end test mass"), 
separated by 2 km or 4 km, form a Fabry Perot resonant op- 
tical cavity with a finesse of approximately 220. Because the 
Michelson interferometer antisymmetric port is held at a dark 
fringe, and because the Fabry-Perot cavities are low-loss, most 
of the light returning from the arms to the beam splitter nom- 
inally exits through the symmetric port of the beam splitter 
back towards the laser. A "power recycling" mirror returns it, 
resonantly, to the interferometer (forming a "power recycling 



cavity"). The average length of the arm cavities is used as a 
frequency reference for the final stage of frequency stabiliza- 
tion [17]. Differential arm cavity length changes result in a 
small amount of light exiting the asymmetric port of the beam 
splitter; this constitutes the gravitational wave signal. The ef- 
fect of the arm cavities and power recycling is to increase the 
sensitivity of the interferometer to gravitational wave signals. 
The arm lengths and arm cavity finesse are optimized to min- 
imize various noise sources. 

The mirrors of the interferometer [20, 2 1 ] are suspended as 
pendulums [22]. Active and passive vibration isolation sys- 
tems [23] are used to isolate them from seismic noise. Various 
feedback control systems are used to keep the multiple opti- 
cal cavities tightly on resonance [24] and well aligned [25], 
and to keep the Michelson interferometer on a dark fringe. 
The LI detector also employed feedforward control to com- 
pensate for microseismic disturbances [26]. When all length 
degrees of freedom are under control and the control systems 
are operating within their linear regime, the interferometer is 
said to be "in lock". During the first few minutes following 
the acquisition of a lock in any individual detector, the in- 
strument typically experiences excess noise due to the ringing 
down of mechanical resonances in the mirror suspensions that 
were excited by impulsive forces applied during the lock ac- 
quisition procedure. After allowing for these resonances to 
damp down, the detector is placed into "science mode"; the 
data collected in science mode are available for gravitational 
wave searches. Science mode continues until the interferom- 
eter loses lock or becomes unstable for any reason. The grav- 
itational wave strain signal (referred to in this paper as the 
gravitational wave data channel) is derived from the error sig- 
nal of the feedback loop used to control the differential length 
of the interferometer arms. A 16 bit analog-to-digital con- 
verter is used to digitize the (uncalibrated) strain signal at a 
rate of 16384 Hz. 

To calibrate the error signal, the response to a known dif- 
ferential arm strain is measured, and the frequency-dependent 
effect of the feedback loop gain is measured and compensated 
for. The laser wavelength and the amplitude of the mirror 
drive signal required to move the interference pattern through 
a fixed number of fringes are used to calibrate the absolute 
scale for strain. The frequency response of the detector is de- 
termined via periodic swept-sine excitations of the end test 
masses. During detector operation, the calibration is tracked 
by injecting continuous, fixed-amplitude sinusoidal excita- 
tions into the end test mass control systems, and monitor- 
ing the amplitude of these signals at the measurement (error) 
point. The calibration procedure, and results, are described in 
more detail in [27, 28]. 



B. The SI run 

By the summer of 2002, all three LIGO detectors were op- 
erating reasonably stably and with reasonable in-lock duty cy- 
cle. As discussed below, the strain sensitivities of all three de- 
tectors were far from their design goals, but were nonetheless 
sensitive to gravitational wave bursts from energetic events in 
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FIG. 1: Typical sensitivities of the three LIGO detectors during the 
SI data run, in terms of equivalent strain noise amplitude density. 
The points are the root-sum-square strain (h ras ) of sine-Gaussian 
bursts for which our TFCLUSTERS analysis pipeline is 50% effi- 
cient, as reported in section V B. 



our Galactic neighborhood. The LIGO Laboratory decided 
that it was an appropriate time for the first Science Run, S 1 . 

The S 1 run consisted of a 408 hour continuous period from 
August 23 through September 9 of 2002, during which data 
were collected from all three LIGO interferometric detectors. 
The state of each of the detectors and the quality of the data 
being logged was continuously monitored through automated 
and manual procedures. As discussed above, in order to be 
sensitive to gravitational waves, the detectors must be in sci- 
ence mode. Environmental disturbances and various instru- 
mental instabilities make it impossible to maintain lock at all 
times, reducing the effective observation time of the run. Dur- 
ing SI, the science mode duty cycles of the three detectors 
were 41.7% for LI, 57.6% for HI, and 73.1% for H2. The 
burst search reported here makes use of the data when all three 
detectors were in science mode simultaneously, comprising 
95.7 hours, or 23.4% duty cycle. 

The strain sensitivity of the LIGO detectors is a strong func- 
tion of gravitational wave frequency. In this analysis, we 
focus on a "detection band" of best strain sensitivity, from 
150 to 3000 Hz. Figure 1 shows amplitude spectra of strain- 
equivalent noise, typical of the three LIGO detectors during 
the SI run. The LIGO design strain sensitivity is also indi- 
cated for comparison. The differences among the three spectra 
reflect differences in the operating parameters and hardware 
implementations of the three instruments; they are in various 
stages of reaching the final design configuration. All detectors 
operated during S 1 at lower effective laser power levels than 
the eventual level of 6 W at the interferometer input. Other 
major differences between the SI state and the final config- 
uration were partially implemented laser frequency and am- 
plitude stabilization systems and partially implemented align- 
ment control systems. Because of these conditions, the strain 
sensitivities of the three detectors were far from the design 
sensitivity (see Figure 1). 



C. Data preselection 



The data processing pipeline described in section III makes 
use of many adjustable parameters that can be tuned to opti- 
mize the search effectiveness. We performed these optimiza- 
tions on a subset of the S 1 data that was reserved exclusively 
for the purpose, and then not used further in the generation of 
scientific results. We called this reserved data set the "play- 
ground" data set. It was chosen to be about 10% of the total 
available triple coincidence data. The choice of which data to 
include was made by hand, to include as much variety of data 
quality as possible. The same playground data set was used 
for both the burst search and the search for inspiralling binary 
neutron stars [12]. This tuning procedure is described in sec- 
tion III. Further, the data processing pipeline analyzed triple- 
coincidence data in six-minute stretches, for convenience in 
data handling. Lock stretches that were less than six minutes 
long, or data in the last < 6 minutes of a longer lock stretch, 
were excluded from further analysis. After exclusion of the 
playground data and these lock stretch boundaries, 80.8 hours 
of triple-coincidence data remain. 

Much effort has gone into improving the stationarity of the 
statistical properties of the detector noise, and understanding 
the noise fluctuations. However, both the detectors' responses, 
and their noise levels, were far from stationary, largely be- 
cause the control systems were not yet completely imple- 
mented. In order to ensure that the data used for this burst 
search are of the highest available quality, we excluded locked 
stretches in which the noise in the gravitational wave chan- 
nel exceeded a pre-determined threshold. The band-limited 
root-mean-square (BLRMS) noise power in the gravitational 
wave channel was monitored continuously in four bands (320 
- 400 Hz, 400 - 600 Hz, 600 - 1600 Hz, and 1600 - 3000 
Hz). Whenever the BLRMS over a six-minute interval for 
any detector in any of these bands exceeded a threshold of 3 
times the 68th percentile level for the entire run (10 times for 
the 320^-00 Hz band), the data from that six-minute period 
were excluded from further analysis. A total of 54.6 hours of 
triple-coincidence data remains after this "BLRMS cut". A 
sufficiently strong gravitational wave burst could trigger the 
BLRMS cut and thereby prevent its own detection; the re- 
quired amplitude is calculated in section V B. 

As discussed in section II A, the response of the detectors 
to gravitational waves was tracked by injecting sinusoidal cali- 
bration excitations into the end test mass control systems. Due 
to technical difficulties, these calibration lines were not reli- 
able or available during some data taking periods. In order 
to ensure that all the data used in this search represent ob- 
servations from detectors with well-understood response, data 
that show no, or anomalously low, calibration lines were ex- 
cluded from further analysis (the "calibration cut"), leaving 
35.5 hours of triple-coincidence data remaining. This is the 
final data sample used to search for gravitational wave bursts. 
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HI. THE DATA PROCESSING PIPELINE 

In the analysis presented here, the purpose of the data pro- 
cessing pipeline is to identify candidate gravitational wave 
events in the data from all three detectors in coincidence. In 
this section, we discuss the procedures and algorithms used 
to identify coincident burst event candidates, the tuning of 
the most important parameters, and the procedures used to 
estimate the accidental coincident burst event rate. The en- 
tire analysis procedure, parameter tuning, event property esti- 
mation, and all other optimizations were developed using the 
playground data (section II C), and frozen before applying the 
analysis to the full SI data set. In the process of analyzing 
the full data set, it became clear that many of the procedures 
and tunings were less than optimal, for a variety of reasons. 
We present the results of this first analysis in this paper, and 
intend to apply improved methods and optimizations (see sec- 
tion VII B) to the analysis of future data sets (which will have 
much greater sensitivity to gravitational wave bursts). 



A. Pipeline overview 

Figure 2 shows, in graphical form, the data processing 
pipeline used in this analysis. Most of the figure is used to 
schematically illustrate various steps in the pipeline of one of 
the interferometric detectors (HI, LI, or H2, generically re- 
ferred to here as IFO-1). The analysis pipelines of the other 
two IFOs (IFO-2 and IFO-3) are not shown in detail because 
they are identical to the first. The first step in the pipeline 
("Band limited RMS & calibration cuts") validates the strain 
channel data used in the analysis; only validated data (section 
II C) taken at times when all three detectors were operating 
simultaneously in science mode are used in this analysis. This 
step establishes the accumulated observation time, or livetime, 
for the analysis. 

The next steps in the pipeline ("Prefiltering & whitening") 
take as input the raw gravitational wave channel data from 
each detector, and prefilter the data stream (section III B). The 
following step ("Burst event trigger generation") searches for 
bursts in the filtered data stream using two different burst de- 
tection algorithms (section IIIC), resulting in a set of event 
triggers at each detector. All data were processed in non- 
overlapping segments that were six minutes long. 

Our pipeline allows for the elimination of event triggers 
that are coincident in time with anomalous events in auxiliary 
channels that monitor the detector and the environment (see 
"Auxiliary channels" path and "Single IFO analysis" gate in 
Fig. 2). The consideration of these potential vetoes will be 
described in section III E. 

Real gravitational wave bursts will cause a nearly simulta- 
neous response in all three detectors, so in the next step we 
require temporal coincidence of single-detector event triggers 
("Multi-IFO analysis" block in Fig. 2). We select as "Event 
candidates" only those combinations of single-detector event 
triggers that are consistent with originating from a single plane 
gravitational wave burst incident on the detector array (section 
III F). 



Event triggers from the single-detector steps of our pipeline 
are mostly due to detector noise. The independence of noise 
triggers at the two LIGO sites is an important assumption in 
this analysis. The largest sources of noise above 100 Hz are 
generated internal to the detectors and are thus uncorrelated. 
Environmental disturbances can produce bursts of noise that 
cause triggers, and these can be coincident between the two 
sites within the ±10 ms gravitational wave travel time if they 
propagate electromagnetically. Auxiliary sensors (e.g., mag- 
netometers) monitor the environment and their inter-site cor- 
relations have been studied. Inter-site disturbances are calcu- 
lated not to be important until the detectors are close to design 
sensitivity, and our studies of S 1 data have not found evidence 
for coincident noise bursts even for H1-H2 where the coin- 
cident location increases the susceptibility to environmental 
disturbances. There does exist some coherence between the 
gravitational wave channels from the three detectors at certain 
frequencies [14], but there is no evidence that this contributes 
to coincident noise bursts. 

Gravitational wave burst events detected at the two LIGO 
sites will be correlated in time. We can evaluate the mean rate 
of background events by measuring the mean rate of events 
that pass our coincident step after we have artificially shifted 
in time all the event triggers identified in one of the detectors, 
for example, L 1 . This background rate estimation is described 
in section IV. 

Finally, to determine the efficiency of the data processing 
pipeline to the detection of strain events incident on the detec- 
tor array we add simulated events, of varying waveform and 
amplitude, to the input data stream and measure the fraction 
identified as event triggers in each detector. Knowing the de- 
tectors' sensitivity to gravitational waves incident from differ- 
ent directions we can combine the results of these simulations 
to determine the mean efficiency for detection of gravitational 
wave burst events incident on the detector array. The effi- 
ciency determination is described in section V. 



B. Prefiltering 

The event trigger generators we employ are designed to 
process data with a white noise spectrum (constant power 
spectral density as a function of frequency). The raw grav- 
itational wave data from all three detectors during SI are 
strongly colored, consisting essentially of randomly fluctuat- 
ing noise with a strongly frequency-dependent power spec- 
trum. These data can be converted to a noise equivalent strain 
signal through a response function which is also strongly 
frequency-dependent, and which is determined through the 
calibration procedure described in section II. The noise also 
contains unwanted features such as spectral lines associated 
with interference from the 60 Hz power mains, mechanical 
resonances in the detector components, and other imperfec- 
tions. 

For the analysis presented here, the data from the gravi- 
tational wave data channel are passed through a linear filter, 
consisting of a 6-th order Butterworth high-pass filter with 
150 Hz cut-off frequency to suppress large noise fluctuations 



7 



Hardware 
injections 



Software 
injections 



LIGO 
Interferometer 
(IFO-1) 



LIGO IFO-3 1 
LIGO IFO-2 I 




Auxiliary channels 

— \ High pass 
— J filtering 



— \ Event-by-event 
veto analysis 



Strain channel 



Band limited RMS 
& calibration cuts 



Veto Triggers 




Prefiltering 
& whitening 



Event Triggers 






Burst event trigger 
generation 



Single IFO 
analysis: 

Veto and event 
Trigger anti 
coincidence 




Multi-IFO 

analysis: 

temporal 
coincidence, 
clustering, 
and frequency cut 



Event 
candidates 



FIG. 2: Schematic outline of the SI bursts analysis pipeline. 



which were apparent at lower frequencies, and a whitening 
filter to flatten the noise spectrum at frequencies above 150 
Hz. Because of the high-pass filtering, we are insensitive to 
Fourier components of a gravitational wave burst below 150 
Hz. The whitening filters are determined using data taken just 
prior to the S 1 run, and are different for each of the three de- 
tectors. No attempt has been made to incorporate the variation 
of the noise power with time, or to otherwise optimize the fil- 
tering. Further, no attempt has been made to remove spectral 
lines from the data stream or suppress their effect on the event 
trigger identification. It is likely that such pre-processing will 
be necessary for future, more sensitive searches with LIGO 
data [29]. 

The impulse response of the prefilter used for this analy- 
sis has a strong ringing, extending to 40 ms. As a result of 
this ringing, the event trigger generation algorithms recognize 
an impulsive event in the strain channel as a cluster of events 
over a long period of time compared to both the sample rate 
and the light travel time between the detectors. This has im- 
portant consequences for the event trigger time resolution and 
the time coincidence of event triggers generated in different 
detectors, as described in section III F below [29]. 



C. Event Trigger Generation 

We use two different techniques to identify event triggers 
from the prefiltered gravitational wave data channel at each 
detector. One technique, which we refer to as SLOPE, is based 
on Refs. [30, 31]. The second technique, which we refer to as 
TFCLUSTERS, is described in Refs. [32,33]. SLOPEandTF- 
CLUSTERS are two different approaches to identifying and 
selecting infrequent transient events that do not share the sta- 



tistical characteristics of detector noise and thus might be of 
gravitational wave origin. These algorithms are implemented 
within the LIGO Data Analysis System (LDAS [34]) environ- 
ment. 

The following discussion of the SLOPE and TFCLUSTERS 
event trigger generators describes and specifies the parameters 
that can be adjusted in order to optimize the performance of 
the algorithms. Some of the parameters can be established 
without reference to the data, since they effectively determine 
the response of the algorithm to the duration (4 ms to 100 
ms) and frequency band (150 to 3000 Hz) characteristics of 
the bursts that are targeted in this search. Others have been 
optimized using the playground data defined in section II C. It 
is assumed that no (or very few) real gravitational wave bursts 
were present in the playground sample. All parameters were 
fixed prior to the processing of the full data set, in order to 
minimize the chance of bias in event trigger generation. 

The parameter optimization, especially the choice of 
thresholds, is guided by competing demands. Lower thresh- 
olds on excess power or amplitude variations result in higher 
rates of event triggers caused by noise fluctuations, but also 
result in higher sensitivity to gravitational wave bursts. The 
criterion we adopted consists of minimizing the upper limit 
for a suite of simulated gravitational wave bursts, described in 
section V. This minimization was applied to the playground 
data set where no triple coincidence event was found after the 
thresholds were fixed. This was consistent with a goal of ob- 
taining a total number of accidental coincident triggers of or- 
der unity, when extrapolated to the remaining 90% of the full 
S 1 data set. 

Nevertheless, the playground data did not adequately repre- 
sent the full S 1 data set, and a variety of additional effects (in- 
cluding the ringing in the prefiltering, as discussed in section 
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III B) resulted in imperfect optimization of the data processing 
pipeline for both event trigger generators [29]. Therefore, the 
resulting number of estimated accidental coincidence events 
was somewhat larger than one, as discussed in section IV C 
and Table I. 



1. SLOPE 

The SLOPE algorithm identifies candidate gravitational 
wave bursts via a threshold on the output of a linear filter 
applied to the prefiltered gravitational wave data in the time 
domain. We choose a filter that is essentially a differentia- 
tor (in time), and trigger on a slope in the data stream which is 
(statistically) inconsistent with expectations from white Gaus- 
sian noise. The SLOPE algorithm is most sensitive when the 
detector noise in the strain channel is whitened. 

The parameters of the SLOPE filter have been tuned so that 
its highest sensitivity is for bursts in which the signal am- 
plitude is increasing linearly with time for ten data samples 
(10 x 61/iscc). The response of the filter to sine waves rises 
with frequency from zero at DC, reaching its first and highest 
maximum at 1 . 1 kHz. Above this frequency, the response of 
the filter falls off, passing through several zeros and secondary 
maxima. Its 3 dB bandwidth is about 1.4 kHz [29]. 

The filter output is searched for extrema indicating the pres- 
ence of bursts. The peak search algorithm compares each suc- 
cessive filter output value with a threshold. If a filter output 
value is found to exceed the threshold, then that point and 
some number of output filter value after the first point exceed- 
ing the threshold are further analyzed. For the analysis con- 
sidered here, 49 output filter values including the point that 
passed threshold are examined, a time interval of 3.0 ms. The 
output filter value having the highest value in this time inter- 
val generates a single trigger. The amplitude of the trigger and 
time of the trigger are written to a trigger database. For this 
analysis, the threshold was fixed and did not adapt to changing 
noise levels [29]. 



2. TFCLUSTERS 

The TFCLUSTERS event trigger generator is a detection 
algorithm which identifies connected regions {clusters) in a 
time-frequency plane where the power is not consistent with 
the expectations for stationary, colored Gaussian noise. The 
TFCLUSTERS algorithm is described in detail in [33], and 
various aspects of its implementation for real data are dis- 
cussed in [32]. The implementation of TFCLUSTERS used 
for our analysis is described below. 

The data from a six minute long segment are first prefiltered 
as described in section IIIB. A time-frequency spectrogram 
is constructed from 2880 periodograms calculated from 125 
ms long non-overlapping subsegments of the six minute long 
segment [29]. 

A first level of threshold is applied to the spectrogram, re- 
sulting in a high-contrast pixelization. 2880 different mea- 
surements of the power are available for every frequency band 



of the spectrogram. Processing one frequency band at a time, 
the power measurements are fit with a Rice distribution [35]. 
Given this fit to the data, the Rice distribution is integrated 
from a power rj to infinity, and rj is varied until the integral 
is equal to a certain pre-defined fraction p. All the pixels of 
the spectrogram with power larger than r\ are then labeled as 
black pixels, while pixels below the threshold are labeled as 
white pixels. The procedure was repeated for all the frequency 
bins in the spectrogram. The number p is called the black 
pixel probability: in the absence of signals, any pixel in the 
spectrogram has, to a good approximation, an equal and inde- 
pendent probability p of being black, in each frequency band. 
Because of this procedure, the effective threshold for black 
pixels varies in response to changing detector noise levels; the 
threshold is "adaptive", as opposed to the fixed threshold em- 
ployed in the SLOPE algorithm. 

The black pixels are then clustered, to look for bursts of 
excess power in a limited region of the time-frequency plane. 
Two levels of clustering are used by TFCLUSTERS, based on 
a study of simulated bursts with varying waveforms. First, a 
cluster is defined as the set of all black pixels which has at 
least one black nearest neighbor (i.e., was touching a black 
pixel by an "edge") in the set. All clusters containing at least 
five pixels are declared significant in this analysis. Second, 
clusters which are not significant according to the latter crite- 
rion are paired together. If the clusters in a pair are closer to 
each other in the time-frequency plane than a certain distance 
threshold, the pair of clusters is declared significant. 

Clusters satisfying the first clustering condition on the raw 
size of a cluster are counted as event triggers. For clusters sat- 
isfying the second clustering condition, generalized clusters 
are formed by linking all the clusters which satisfy the dis- 
tance thresholds, and these generalized clusters are counted as 
event triggers. For each event trigger, the time and frequency 
intervals over which the cluster extends, the total amount of 
power in the cluster, and the number of pixels it contains, are 
stored in a database. The total power in each cluster is a mea- 
sure of the signal-to-noise ratio for the burst event. It is cal- 
culated without reference to the response of the detector to 
gravational wave bursts, so its relationship to the strength of 
the burst depends on the detector and frequency band. 

The black pixel probability p is tuned as described above. 
The values are different for the three different detectors and 
vary from 0.02 to 0.05. The total power in the cluster 
is required to exceed a pre-determined threshold in post- 
processing; this is effectively a cut on the signal-to-noise ratio 
for the burst event. The threshold on the power is the same 
for all three detectors, in order to obtain rates for false (noise) 
triggers which are roughly the same for all three detectors. 



D. SLOPE and TFCLUSTERS event triggers 

Figs. 3 and 4 show histograms of SLOPE and TFCLUS- 
TERS event triggers before and after the application of the 
BLRMS and calibration cuts described in section II C. The 
horizontal axis in these histograms is a measure of the ampli- 
tude or power of the excess signal identified by the SLOPE 
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FIG. 3: Histogram of SLOPE event triggers from the three LIGO 
detectors, before and after the BLRMS and calibration cuts. 



FIG. 4: Histogram of TFCLUSTERS event triggers from the three 
LIGO detectors, before and after the BLRMS and calibration cuts. 



or TFCLUSTERS algorithms, respectively. These measures 
are indeed proportional to the true amplitude or power of a 
detected gravitational wave burst, as demonstrated in section 
V. However, no information about the detectors' calibrated 
response functions is used in forming these measures, so the 
proportionality constant is different for different waveforms, 
detectors, and data epochs (and is taken into account in the 
evaluation of the detection efficiency, section V). The lower 
limits on the horizontal axis in these histograms correspond 
to the threshold applied to that event trigger for input into the 
next step in the data processing pipeline (triple coincidence). 



E. Auxiliary Channel Vetoes 

Environmental disturbances and detector instabilities could 
also produce event triggers. We collect data in a large number 
of auxiliary channels which monitor the detector and the envi- 
ronment, in order to look for time-coincident bursts and thus 
form vetoes for such false triggers. Our pipeline has the capa- 
bility to search for such bursts in auxiliary channels, and veto 
an event trigger if it is time-coincident with such a burst. Engi- 
neering runs performed prior to the S 1 run indicated that such 
vetoes could be very efficacious, reducing the rate of false 
event triggers with minimal loss of livetime, due to clearly 
identifiable instabilities in the detectors. However, once these 
instabilities were identified, they were eliminated through im- 
proved instrumentation, resulting in much improved stability 
during S 1 . After careful study, no vetoing criteria using aux- 
iliary interferometer and physical environment monitor chan- 
nels are found to be especially efficacious in the S 1 data, for 



this burst search. The most promising vetoing channels in the 
SI data are interferometer sensors that are closely related to 
the gravitational wave channel. While we investigated a num- 
ber of such channels and methods for identifying veto criteria, 
in the end we concluded that further study was needed before 
any of these could safely be used to exclude data from analy- 
sis. Further, employing the identified vetoes would have made 
a negligible difference in the results of this analysis. Thus, in 
this analysis, we apply no vetoes based on auxiliary channels. 

F. Coincidence 

The final stage of our data processing pipeline brings to- 
gether the event triggers generated by a particular event trig- 
ger generator (either SLOPE or TFCLUSTERS) and assem- 
bles a smaller list of coincident event trigger triplets. Each 
triplet consists of an event trigger from each detector that oc- 
cur within an interval consistent with their origin in a single 
gravitational wave burst. These triplets are the event candi- 
dates that form the basis for our determination of bounds on 
the rate of gravitational wave bursts incident on the Earth. 

Temporal coincidence is the most obvious application of 
coincidence for selection of gravitational wave events and ex- 
clusion of noise events. The LIGO detectors are approxi- 
mately co-aligned and co-planar. As a result, they all sense 
approximately the same polarization of any incident gravita- 
tional wave. Correspondingly, all estimated parameters of the 
burst (such as strain amplitude and frequency band) should be, 
up to uncertainties in the estimation, the same for all three de- 
tectors (after accounting for the differences in the detectors' 
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sensitivities). In the analysis presented here we require tem- 
poral coincidence (to an appropriate precision) for both the 
SLOPE and the TFCLUSTERS pipelines. Additionally, TF- 
CLUSTERS events are also characterized by frequency infor- 
mation; we require consistency between the frequency bands 
in a coincident triplet (section III F 3). No attempt is made to 
require coincident event triggers to have consistent amplitudes 
and waveforms [29]. 

In the remainder of this section we describe in greater detail 
the elements of the data processing pipeline coincidence step. 



1. Temporal coincidence 

Gravitational waves arrive at the Earth as plane waves. 
Since gravitational waves are assumed to propagate at the 
speed of light, the interval between event triggers in the differ- 
ent detectors should be no more than the greater of the light 
propagation time between the detectors and the uncertainty 
in the arrival time determination of a prototypical burst asso- 
ciated with the event trigger generator. Different timing un- 
certainties are associated with different event trigger genera- 
tors. Correspondingly, we use different window durations for 
SLOPE and TFCLUSTERS. Given a window, we compare the 
start times of the event triggers generated in each of the three 
detectors. We form an event trigger triplet, or triplet for short, 
from all combinations of HI, H2 and LI events whose start 
times all lie within the window duration. 

As described in section III B the input to the event trigger 
generators is processed through a high-pass filter that rang 
strongly. As a result of this ringing, impulsive events lead 
to a train of multiple SLOPE triggers, with a total duration of 
approximately 40 ms. We add 10 ms to this, corresponding to 
the light travel time between detectors, to determine a 50 ms 
window for temporal coincidence of SLOPE events [29]. 

As described in section IIIC2, TFCLUSTERS was tuned 
to a natural time resolution of 125 ms, much larger than the 
light travel time between the detectors. On the basis of studies 
which indicated a larger range of trigger time differences for 
simulated signals, we expanded this and use a 500 ms window 
to determine triplets of temporally coincident TFCLUSTERS 
events [29]. 



where "something unusual" has happened simultaneously at 
multiple detectors, we want to cluster these sets of closely 
spaced events. 

Clustering takes place only after the time coincidence step. 
We require a minimum separation in time between distinct 
coincident trigger triplets, of 0.5 seconds; triplets that are sep- 
arated in time by less than this amount are clustered together 
into one clustered event triplet (event candidate). The choice 
of the clustering window is based on the study of noise trig- 
gers and simulated bursts (section V B). In the TFCLUSTERS 
pipeline, 0.5 seconds is the width of the coincidence win- 
dow between triggers from the three detectors. In the SLOPE 
pipeline, the coincidence window of 50 ms is too small a sep- 
aration to avoid ambiguities in the definition of clusters and 
in the event counting, so we use 0.5 seconds for consistency 
with the TFCLUSTERS pipeline [29]. 

All triggers in the cluster are assumed to originate from 
one burst event. Guided by simulation studies (section V B), 
the start time, frequency band, and amplitude or power of the 
event is taken to be that of the trigger with the largest ampli- 
tude or power in the cluster. 



3. TFCLUSTERS frequency cut 

For TFLUSTERS we apply one more criterion in the coin- 
cidence step of the pipeline. A triplet of event triggers that 
arises from a single gravitational wave burst incident on all 
the detectors should have consistent values for the estimated 
parameters of the burst. TFCLUSTERS characterizes each 
burst event trigger by its bandwidth: the low and high fre- 
quency bound (/i ow , /high) of the cluster identified in the time- 
frequency plane. When multiple triggers from one detector 
are clustered in time as described in section IIIF2 above, 
the inclusive frequency band for that clustered event trigger 
is formed. For TFCLUSTERS triggers only, we require that 
the frequency bands of the clustered event triggers from each 
detector in the triplet either overlap, or are separated in fre- 
quency space by no more than a fixed window of A/ = 80 Hz, 
based on studies of the simulations described in section V B. 



2. Clustering 

The next step in the multiple-detector coincidence analysis 
is to cluster the events from each detector (this is unrelated to 
the pixel clustering that forms the heart of the TFCLUSTERS 
event trigger generation, section IIIC2). Both the TFCLUS- 
TERS and the SLOPE event trigger generators often associate 
several event triggers with the same "burst" feature. For in- 
stance, the ringing of a 1 ms Gaussian due to the detector 
response and the prefiltering of the data (section IIIB) can 
produce several closely spaced event triggers. TFCLUSTERS 
often associates multiple triggers with the same broadband 
event, all with the same start time but different frequency. 
Since we are interested in the identification of time intervals 



IV. BACKGROUND AND SIGNAL RATES 

The data processing pipeline (section III) generates back- 
ground event triggers originating in noise level fluctuations in 
the detectors, due to random processes or environmental or 
instrumental disturbances. Our primary means to reject such 
background event triggers is temporal coincidence between 
the three detectors in the LIGO array (section III F). To the 
extent that noise fluctuations in each of the detectors are ran- 
dom, uncorrelated, and follow Poisson statistics, the primary 
background comes from accidental coincident events, and the 
accidental triple-coincidence rate can be predicted from the 
observed instantaneous single-detector event rates. 



11 



A. Background estimation 

We have chosen to tune our event trigger generators (using 
the playground data sample) so as to produce an estimated 
accidental triple-coincidence rate of one event over the entire 
S 1 observation time, as discussed in section III C. 

Again assuming no correlations between noise fluctuations 
in the three detectors, we can indirectly measure the rate 
of accidental triple-coincident events from triple-coincidence 
rates when artificial time shifts are introduced between single- 
detector event triggers. 

Such time shifted triple-coincidence events are free of con- 
tamination from true gravitational wave bursts (assuming that 
such bursts are rare), and thus are an unbiassed estimate of 
the accidental triple-coincidence rate. The distribution in the 
number of time shifted triple-coincidence events should fol- 
low a Poisson distribution. These distributions can be fitted to 
obtain the expected number of background events for use in 
our statistical analysis. 

The time shifts should be larger than the maximum dura- 
tion of a real (noise-induced or gravitational wave-induced) 
detectable burst, or else the events will be correlated and 
will not obey Poisson statistics. The time shifts should also 
be shorter than the typical time scale over which the single- 
detector event rates vary substantially, so that the number of 
events for different time shifts will be Poisson distributed for 
a quasi-stationary process. 

To establish a lower limit on the time shift required to 
ensure uncorrected noise event triggers, we histogram the 
time delay between consecutive events in the three detectors, 
shown in Fig. 5 for the TFCLUSTERS event trigger genera- 
tor (the distributions are similar for the SLOPE event trigger 
generator). The distributions of delay times follows the ex- 
pected exponential form for delay times exceeding 8 seconds 
(vertical dashed lines), for all three detectors. Any residual 
auto-correlations present in the data will rapidly decay for de- 
lay times exceeding 8 seconds, and in the case of many (TV) 
time-shift experiments, their potentially biased contribution to 
the Poisson estimate reduce as 1/N, 

The assumption that noise fluctuations are uncorrected 
between detectors is questionable for the two detectors co- 
located at the Hanford site [36], HI and H2. Indeed, there ex- 
ists evidence for short-term, narrow-band correlations in the 
noise power between the HI and H2 detectors associated with 
power line harmonics, as well as correlations between LI and 
HI or H2 associated with harmonics of the data acquisition 
buffer rate [14]. The power line harmonics integrate away 
over long time scales, and the data acquisition buffer rate har- 
monics only appear after long integration times. It is the short 
term correlations that concern us here. We have found no de- 
tectable evidence of short term correlated noise fluctuations 
associated with these sources of narrow-band correlations. 

In order to account for any potential correlations in noise 
fluctuations between HI and H2, we have performed our time- 
shifted coincidence measurements by shifting the time be- 
tween event triggers found in the LI data and those found in 
the HI and H2 data, while keeping zero time shift between HI 
and H2. 
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FIG. 5: Histograms of the time delay between consecutive events 
in the TFCLUSTERS event trigger generation, for the LI, HI, and 
H2 detectors. The curves are components of fits to the distributions 
that incorporate the expectations for short time delay correlations and 
long time delay random, uncorrected events. The vertical dashed 
lines indicate the time delay beyond which consecutive events are 
consistent with being uncorrelated. 



We have performed multiple time-shift experiments with 
the SLOPE and TFCLUSTERS event trigger generators be- 
tween the Livingston and Hanford sites. The resulting number 
of time-shifted triple coincident events from 24 such experi- 
ments in the [—100, 100] second range with 8 second steps 
are shown in Figs. 6 and 7 for TFCLUSTERS and SLOPE, re- 
spectively. The distributions of background events for the 24 
non-zero time shifts (lower plots in Figs. 6 and 7) are fitted 
with Poisson predictions and are found to be consistent with 
the expectation from Poisson statistics. Averages and Poisson 
mean values for different step and window sizes vary by less 
than 0.5 events. 

In estimating the background rate, we have considered time 
shifts between 8 seconds (to avoid correlated events; see 
Fig. 5) and 100 seconds (to minimize dependence on any 
non-stationarity in the instantaneous event rate). These time- 
shift measurements yield estimates of the number of acci- 
dental triple-coincidence (background) events /i b for the TF- 
CLUSTERS and SLOPE event triggers. Because these mea- 
surements can be made with many, statistically-independent 
time shifts, the resulting statistical errors on these estimates 
are small. 

The presence of any remaining non-stationarity in the back- 
ground event rate, however, will result in errors in the back- 
ground rate estimation. In fact, the instantaneous event trigger 
rate is observed to vary for both event trigger generators. The 
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within that segment, the singles rate varies markedly, making 
it difficult to reliably estimate the background rate. For this 
reason, we choose not to use the SLOPE pipeline to set a limit 
on the rate of gravitational wave bursts [29]. 

It should be noted that before these coincidences and back- 
ground rates were available, we decided to set our upper lim- 
its using the results from the event trigger generator which 
yielded the better efficiency for detecting gravitational wave 
bursts, as measured by our simulations. For almost all wave- 
forms this turned out to be the TFCLUSTERS pipeline (see 
section V D). Thus, even if the background rate (and thus the 
rate of excess triggers) from the SLOPE pipeline could be re- 
liably estimated, the primary results from this search would 
still be based on the TFCLUSTERS pipeline. 

B. Signal candidate estimation 



FIG. 6: Time-shifted triple coincident events from TFCLUSTERS 
event triggers, as a function of an artificial time shift introduced be- 
tween the Hanford (LHO) and Livingston (LLO) sites. Top: Num- 
ber of events versus time shift, in 8 second steps; the point at zero 
time shift is the number of true triple coincident events. Bottom: 
Histogram of the number of time-shifted coincident events, with the 
Poisson fit overlaid (the zero time shift point is excluded). In both 
plots, the error bars are Poissonian. 
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FIG. 7: Same as Fig. 6, but for the SLOPE event trigger generator. 

variability of the trigger rate is sufficiently severe for SLOPE 
that the background rate estimation is deemed unreliable. Be- 
cause of the fixed (non-adaptive) threshold employed in the 
SLOPE algorithm, the trigger rate of the individual interfer- 
ometers varies by more than a factor of one thousand, some- 
times on timescales of 10 seconds or less. All five events of 
the zero time lag coincidences and most of the time-shifted 
coincidences in Fig. 7 come from a single 360-second seg- 
ment, corresponding to a coincidence of data segments from 
three interferometers with very high burst singles rates. Even 



An excess in the number of coincident (zero-time-shift) 
events over the estimated background can be estimated statis- 
tically. Here we make use of the unified approach of Feldman 
and Cousins [37]. This approach provides an unambiguous 
prescription for establishing a statistical excess of signal can- 
didate events at a specified confidence level (that is, a lower 
limit to the confidence interval that is greater than zero). How- 
ever, as discussed in section I, we have not yet characterized 
our detectors and data analysis procedures sufficiently well to 
claim that any such excess is a detection of gravitational wave 
bursts. We therefore use only the upper endpoint on the confi- 
dence interval for the number of signal candidate events to set 
an upper limit on the rate of gravitational wave bursts. 

Starting from an observed number of events n and an esti- 
mate of the number of background events hb, we build confi- 
dence bands for the number of signal events [is according to 
the formula: 

P( n o)ns+»B > at, (4.1) 

where p(n)^ s +^ lB is the Poisson probability density function 

PWw+w = (Ms + Ms)" : ■ (4.2) 

n! 

The sum extremes, n\ and n<i, are chosen according to a like- 
lihood ranking principle [37]. In our implementation, we as- 
sume both signal and background are Poisson distributed. We 
report confidence bands for a = 90%, 95%, and 99%. 

We account for the statistical error on the background es- 
timation following the method described in [38, 39], where a 
Gaussian background uncertainty is folded in the formulation 
of the probability density function. We replace p(n) MS _|_^ B in 
Eqn. 4.1 with 

qHhs+hb = -/== — / P( n )ns+n' B e dn' B 
V2tt(Tb Jo 

(4.3) 

where erg is the estimated background error. This marginal- 
ization, performed through a Monte Carlo calculation, is used 
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in the construction of confidence bands for the estimated back- 
ground /ig ± erg. 

C. Event rate bound 

Table I shows the number of coincident events, the esti- 
mated number of accidental coincident events (background), 
and the confidence bands that we find at the 90%, 95%, and 
99% confidence levels on the number of excess events and 
the event rate (over 35.5 hours of SI observation time), using 
the TFCLUSTERS event trigger generator pipeline. The up- 
per bounds of the confidence bands are taken to be the upper 
limit on the number of signal events, at that confidence level. 
At the 90% confidence level, the search yields an upper limit 
of 2.3 events in 35.5 hours. As discussed in section IV A, be- 
cause of the variability of the event trigger rate in the SLOPE 
pipeline, we choose not to use it to set a limit on the rate of 
gravitational wave bursts. 

TABLE I: Confidence bands on the number of excess events in the 
SI run (35.5 hours of observation time) from the TFCLUSTERS 
pipeline. 

Coincident events 6 
Background 10.1 ±0.6 

90% confidence band 0-2.3 
95% confidence band — 3.5 
99% confidence band — 5.9 



Given the estimated backgrounds from the time shift anal- 
yses, the number of TFCLUSTERS events at zero time lag is 
somewhat low (Table I and Fig. 6). None of the events de- 
tected by SLOPE were detected by TFCLUSTERS. This is not 
in itself surprising, since the two event trigger generators have 
different sensitivities to different waveforms, but it does indi- 
cate that none of the events were far above threshold for that 
trigger generator, since the largest differences in efficiency be- 
tween the two event trigger generators was approximately a 
factor of 6 (section V D below). The probability of obtaining 
six or fewer TFCLUSTERS events, given our estimated back- 
ground, is approximately 12%. We found no reason to sus- 
pect any systematic errors in our background estimate for this 
pipeline. Alternative methods of estimating the background 
(simple estimates based on the average singles rates and the 
coincidence window, time shift analyses where all three de- 
tectors are shifted as opposed to holding H1-H2 fixed at zero 
delay) did not give significantly different background rates. 

V. EFFICIENCY DETERMINATION 

In order to interpret our bound on the observed rate for co- 
incident gravitational wave bursts, we study the response of 
the LIGO detectors and our analysis pipeline to simulated sig- 
nals with varying waveforms, durations, bandwidth, and peak 
amplitudes. The simulated signals were injected into the grav- 
itational wave data stream from each of the three detectors, as 



far upstream in the pipeline as was practical (after data acqui- 
sition and ingestion into the LIGO Data Analysis System). 

The same data that were used to search for coincident bursts 
(section II C) were also used for these simulations; and, for 
the purposes of these simulations, these data are assumed to 
consist entirely of noise (no real gravitational wave bursts 
present). Approximately 20% of the SI data was used for 
these simulations, spanning the entire data run uniformly, in 
order to fairly represent the noise and detector sensitivities 
throughout the run. 

We present the results for the efficiency determinations for 
both TFCLUSTERS and SLOPE event trigger generators, even 
though (as noted in section IV A) we do not use SLOPE to 
derive a final limit on the rate of gravitational wave bursts. 

A. Waveforms 

The astrophysical origin, and waveform morphology, of the 
gravitational wave bursts we search for in this work are a pri- 
ori unknown. A broad range of signal waveforms were con- 
sidered. These include astrophysically-motivated waveforms, 
such as the results of supernova simulations [40, 41], as well 
as ad-hoc waveforms such as Gaussians, damped sinusoids, 
sine-Gaussians, Hermite-Gaussians, and others. Guided by 
the simulations in [40, 41], we have endeavored to be sensi- 
tive to any waveform that adds excess power (over that of the 
detector noise) in the LIGO SI sensitivity band (150 to 3000 
Hz), with durations between 4 ms and 100 ms. 

In order to evaluate our sensitivity to such bursts, we must 
model the waveforms in some general way. For the results 
presented here, we have chosen to focus on two classes of ad 
hoc waveforms, which we regard as "surrogates" for real as- 
trophysical signals. The first are broad-band, limited-duration 
Gaussians of the form 

h(t + t ) = /i cxp(-£ 2 /V 2 ) , (5.1) 

with varying peak amplitude ho, peak time to, and duration 
t (Fig. 8). The second are narrower-band, limited-duration 
sine-Gaussians of the form 

h(t + t ) = h sin (27r/ i) cxp (-t 2 /t 2 ) . (5.2) 

The duration of the sine-Gaussians were chosen to be t = 
2/io- Their Fourier transforms h(f) span a (Gaussian) fre- 
quency band of Of = fo/Q centered about the central fre- 
quency /o, where Q = fo/<Jf = V^ttt/o ~ 8.9. We 
have chosen eight different central frequencies, spaced log- 
arithmically, and spanning the LIGO sensitivity band: /o 
100. 153, 235, 361, 554, 850, 1304, 2000] Hz (Fig. 9). 
Our analysis pipeline detection efficiency depends on the 
burst duration, frequency band, and some measure of burst 
"strength"; it does not depend strongly on the precise wave- 
form morphology. In order to facilitate comparison of the 
burst strength with the detectors' equivalent strain noise, and 
with burst waveforms having similar properties but different 
detailed morphologies, we define two useful measures of the 
burst strength. The root-sum-square (rss) amplitude spectral 
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FIG. 8: Gaussian waveforms, with varying duration, as described in 
section V A. Top: Frequency spectrum. Bottom: Time series. 
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FIG. 9: Sine-Gaussian waveforms, with central frequency fo evenly 
spaced in log(f), and Q ~ 9, as described in section VA. Top: 
Frequency spectrum. Bottom: Time series for a sine-Gaussian with 
fo = 554 Hz. 



density for such bursts, in units of dimensionless strain per 
root Hz, is defined by: 




7r/2r/io (Gaussians); 



(5.3) 



(5.4) 



<5/(4v / 7r/o)^o (sine-Gaussians). (5.5) 



The characteristic strain amplitude, in units of dimensionless 
strain, is defined by [42]: 

Khar = fc\h(fc)\; (5-6) 

= Vn(fcT) exp (-(tt/ c t) 2 ) h (Gaussians);(5.7) 



= Vn(f c T/2)h (sine-Gaussians, f c = f ). (5.8) 
Here, h(f) is the Fourier transform of h(t), defined by 

/>00 

h(f) = / h(t) e - l2wft dt, (5.9) 



and f c is a characteristic frequency (typically, either the fre- 
quency at which h(f) peaks, or the frequency where the de- 
tector is most sensitive). For our sine-Gaussians, we choose 
f c to be the central frequency fo = 2/r; for Gaussians, we 
choose f c to be the frequency at which all three LIGO detec- 
tors had approximately best sensitivity during SI, f c ~ 300 
Hz (see Fig. 1). 



B. Simulations 

In order to add the simulated signal (in units of dimension- 
less strain) to the raw detector data (in units of ADC counts), 
we must convert, or filter, the signal using the detector re- 
sponse function (in counts per strain) obtained through the 
calibration procedure described in section II. The simulated 
signals, padded with zeros to minimize edge effects, are fil- 
tered through the detector response function in the Fourier do- 
main, yielding a time series in ADC counts that can be added 
directly to the raw gravitational wave data stream at the begin- 
ning of the data processing pipeline. These simulated signals 
can be injected at any chosen point in time, and with any cho- 
sen amplitude. The uncertainty in the calibration information 
is the largest source of systematic error in this analysis (sec- 
tion VIC). 

For each waveform, we evaluate the efficiency for detec- 
tion through each of the three LIGO detectors and analysis 
pipelines, as a function of h rss (defined in equation 5.3), as- 
suming optimal wave direction and polarization. Approxi- 
mately 80 simulations are performed for each combination of 
waveform, h rss , detector, and event trigger generator, using 
data spanning the SI run. In Figs. 10 and 11 we plot detec- 
tion efficiencies and average signal strengths for the TFCLUS- 
TERS and SLOPE event trigger generators, respectively. Al- 
though our event trigger generators do not necessarily trig- 
ger on excess power, we find that the "strength" of the signal 
reported by either event trigger generator (the SLOPE "am- 
plitude" or the square root of the TFCLUSTERS "power") is 
proportional to the actual amplitude of the injected signal over 
a broad range of h rss . This is illustrated by the lower plots in 
Figs. 10 and 1 1 (for one particular waveform). We emphasize 
that the "power" or "amplitude" that is plotted in Figs. 10 and 
11, respectively (and in Figs. 4 and 3), are purely algorithm- 
dependent quantities which are compared with thresholds to 
define event triggers; they are not designed to be true measures 
of the burst power or amplitude, and they will be different for 
different detectors and waveforms. 

As expected, the efficiencies are essentially 100% for large 
values of h rss , consistent with noise and thus 0% efficiency 
for small h rss , and transitioning smoothly over a narrow in- 
termediate range of h rss . The time window used to associate 
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a TFCLUSTERS event trigger around the time of the injec- 
tion (0.5 s) is larger than for a SLOPE event trigger, so the 
observed efficiency for small h rss waveforms appears larger 
in Fig. 10 than in Fig. 11. The results in both cases are em- 
pirically found to be well fitted to simple sigmoid curves in 
log w (h rss y. 



s{hrss) 



i 



1 + e~( lo Sio hrss~b)/a 



(5.10) 



where b = log 10 h rss i n determines the strain per root Hz 
h rss 1/2 at which the efficiency is equal to 1/2, and a governs 
the width of the transition from to 1 in log 10 (h rss ). It is 
specific to a given waveform, detector, event trigger generator, 
and data epoch. All fits resulted in good fit quality, except at 
the smallest values of h rss , where noise triggers dominate; we 
exclude such triggers from our definition of "efficiency", and 
use the sigmoid fits to extrapolate to zero efficiency at small 
values of h rss . Examples of sigmoid fits are shown in Fig. 10 
and Fig. 11. 




h rss [ Hz 12 ] 



FIG. 10: The response of the TFCLUSTERS event trigger generator 
to Gaussian bursts with r = 1 ms, embedded in S 1 data, as a function 
of the root-sum-square strain h rS s- Upper plot: Average burst detec- 
tion efficiency. The efficiencies were evaluated through simulations 
of burst waveforms with optimal wave direction and polarization, in- 
jected into SI data. The simulated data points are fitted to sigmoid 
curves, shown, in the region where the efficiency is not dominated 
by random noise triggers. The curve for the triple-coincidence is the 
product of the single-detector efficiency curves, and can be directly 
compared with the triple-coincidence simulation data points. Lower 
plot: Average detected signal strength for each of the three LIGO 
detectors. 

The efficiencies shown in Figs. 10 and 11 do not remain 
at unity to arbitrarily high amplitudes. A sufficiently strong 
gravitational wave could trigger the BLRMS cut and in that 
way effectively prevent its own detection in this search. The 
most susceptible band for such a possibility is the 320-400 




h rss [ HZ" 1 ' 2 ] 



FIG. 11: This is the same as Fig. 10 but obtained with the SLOPE 
event trigger generator. 



Hz band in LI, where, for example, a 361 Hz sine-Gaussian 
with h rss ^6 x 10~ 18 could begin to trigger the BLRMS cut. 
Signals centered at other frequencies or those with broader 
bandwidths would require a higher amplitude. We estimate 
that a loud supernova [40, 41] at 7 pc, or a 3+3 solar mass bi- 
nary neutron star inspiral at 300 pc could begin to trigger the 
BLRMS cut. (Note however that the well-defined waveform 
of the latter makes a template -based search [12] a more sensi- 
tive method for detecting such waves.) The need for this cut in 
the data was driven by the nonstationarity of the noise in the 
detectors during S2, and detector improvements are expected 
to reduce our use of such cuts in the future. 

In Fig. 1 we compare the value of h rss for which our sim- 
ulations of sine-Gaussian waveforms (at optimal wave direc- 
tion and polarization) yield 50% efficiency (averaged over the 
entire SI run), shown as circles, with the detectors' (typi- 
cal) equivalent strain noise. These 50% efficiency points are 
roughly an order of magnitude larger than the equivalent strain 
noise. 



C. Average over direction and polarization 

The response of a LIGO detector to an incident gravita- 
tional wave burst depends on the wave direction and wave 
polarization relative to the detector axes, and is referred to 
as the detector's antenna pattern [15]. The only effects of the 
wave direction and polarization are to modify the amplitude 
of the detected wave and the relative arrival times at the de- 
tectors. Since we have evaluated the detection efficiency for 
each detector as a function of the root-sum-square strain of the 
wave at optimal direction (directly overhead) and polarization 
(aligned with the detector axes), it is straightforward to eval- 
uate the efficiency at arbitrary direction and polarization. We 
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choose to consider a population of sources distributed isotrop- 
ically in the sky, with random linear polarization. We thus 
evaluate the detection efficiency averaged over direction and 
polarization, as a function of intrinsic strain per root Hz inci- 
dent on the Earth. 

In order to evaluate the efficiency for coincident detection 
by non-colocated detectors, we assume that the detection ef- 
ficiency is a measure of a random process, uncorrected be- 
tween detectors. Further, the difference in arrival times at the 
different detectors is small compared to the time coincidence 
window employed (section III F 1). Therefore, the efficiency 
for triple coincidence can be expressed as the product of effi- 
ciencies for the three LIGO detectors evaluated at the appro- 
priate peak amplitude for each. Under this assumption, the ef- 
ficiency for coincident detection by all three LIGO detectors, 
averaged over wave direction and polarization, can be evalu- 
ated numerically. The results of this procedure are shown in 
Figs. 12, 13, and 14. 
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FIG. 12: Burst detection efficiency for Gaussian bursts with r = 1 
ms, as a function of h ras , for each of the three LIGO detectors, and 
for the triple coincidence, using the TFCLUSTERS event trigger 
generator. The lighter grey curves are the same as the curves in 
Fig. 10. The darker curves to the right of them are the result of av- 
eraging the efficiency curves over wave directions and polarizations 
(denoted by (■ ■ ■) in the legend) as described in section VC. 
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FIG. 13: Burst detection efficiency for triple coincidence as a func- 
tion of h r3S , using the TFCLUSTERS event trigger generator, aver- 
aging over wave directions and polarizations, for six different wave- 
forms: GA refers to the Gaussians defined in Eqn. 5.1 and SG to the 
sine-Gaussians defined in Eqn. 5.2. 
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FIG. 14: Same as Fig. 13, for the SLOPE event trigger generator. 



The single-detector efficiencies will be independent only if 
there are no significant noise correlations between the detec- 
tors. We have compared the results for coincident detection to 
direct simulations of triple coincidence through the full three- 
detector analysis pipeline (under the simplifying assumption 
of neglecting the difference in antenna pattern response be- 
tween the two sites), and found good agreement (see triple 
coincidence data points in Figs. 10 and 11); no evidence for 
burst detection efficiency correlations between the detectors 
has been observed. 

A crucial test of the accuracy of this simulation proce- 
dure is the comparison of signals injected into the data stream 
with software, with signals injected directly into the end test 
masses of the interferometers. The comparison of these "hard- 



ware" burst injections with the (much more numerous) soft- 
ware injections provide a test of the detector response, the cal- 
ibration information, the data acquisition, and the entire soft- 
ware analysis chain, including the software simulations used 
to evaluate the efficiency, as described above. 

Hardware injections of sine-Gaussians with varying /o and 
h rss were performed during brief periods, just prior to the be- 
ginning and just after the end of the S 1 run. Due to time con- 
straints, only a limited number of hardware injections were 
performed. As a result, the comparison with software injec- 
tions is somewhat qualitative. The detection of these hardware 
injections through the analysis pipeline was found to be con- 
sistent with expectations from the software injections. 
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D. Sensitivity to modeled bursts 



We can use the efficiency determined by simulated signal 
injections, discussed in sections V B and V C, to estimate the 
weakest signal we could have seen in the search described in 
this article. The efficiencies of each of our two event trigger 
generators for several different waveforms is shown in Figure 
13 and Figure 14. The sensitivity at 50% efficiency, for a 
variety of Gaussian and sine-Gaussian waveforms, is shown 
in Table II in terms of h rss , and in Table III in terms of h c har- 



TABLE II: Sensitivity to various waveforms in the SI run from TF- 
CLUSTERS and SLOPE pipelines for triple-coincidence detection, 
averaged over source direction and polarization. The sensitivity is 
given in terms of h rss (Eqn. 5.3, units of Hz" 1 / 2 ), at 50% efficiency 
(h rss i/2)- A 20% uncertainty associated with calibration (section 
VI C) is not included. 
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TABLE III: Sensitivity to various waveforms in the S 1 run from TF- 
CLUSTERS and SLOPE pipelines for triple-coincidence detection, 
averaged over source direction and polarization. The sensitivity is 
given in terms of h h ar (Eqn. 5.6, dimensionless strain), at 50% effi- 
ciency. A 20% uncertainty associated with calibration (section VI C) 
is not included. 
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VI. INTERPRETED RESULTS 

A. Exclusion in rate versus strength plane 

The results of our search can be used to set limits on mod- 
els of ensembles of gravitational waves arriving at the earth. 
Fig. 15 shows the upper limits that we set, using the TFCLUS- 
TERS event trigger generator, as expressed in the plane of 
event rate versus h rss . The top figure is for the case of 1 ms 
and 2.5 ms Gaussian bursts, and the lower figure is for sine- 
Gaussian bursts with central frequency of 361, 554, 850 and 
1304 Hz. 

As discussed in sections I and V, these limits are given in 
terms of an ensemble of waves of equal amplitude, incident 
on the earth from all directions and with all (linear) polariza- 
tions. This ensemble is not motivated by astrophysical consid- 
erations, but is nevertheless useful in characterizing the per- 
formance of the search, and it can be compared with similar 
limits obtained by resonant bar detector collaborations [1,2]. 




FIG. 15: Rate versus h r sa for detection of specific waveforms us- 
ing the TFCLUSTERS event trigger generator. The region above 
and to the right of the curves is excluded at 90% confidence level or 
greater. The effect of the 20% uncertainty in the detector response 
is included. Top: For Gaussians with r = 1.0 ms and r = 2.5 ms. 
Bottom: For sine-Gaussians with Q — 9 and central frequency fo = 
361,554, 850 and 1304 Hz. 

The curves in Fig. 15 are constructed by dividing our ob- 
served event rate upper limit by the efficiency curves shown 
in Fig. 13. In the limit of large wave amplitudes h rss where 
our efficiency is essentially unity, the upper limit is indepen- 
dent of amplitude, at a level given by the burst rate limit pre- 
sented in section IV C. The limit at all amplitudes with lower 
efficiency is given by that burst rate limit, multiplied by the 
inverse of the efficiency. 

The meaning of the upper limit curve can be understood by 
imagining that one is interested in the flux of 1 ms Gaussian 
gravitational wave bursts at a particular amplitude. Reading 
the value of the curve at that amplitude gives the 90% con- 
fidence level upper limit on the rate of such bursts with h rss 
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equal to or greater than that amplitude. Similarly, one can 
use these curves to determine the 90% confidence level up- 
per limit on the amplitude of bursts which are incident at a 
specified rate. 



B. Comparing results from the two pipelines 

As discussed in section IV A, the variability in the SLOPE 
background led us to choose not to use it to set an upper 
limit on the gravitational wave burst rate. However, the two 
event trigger generators use very different and complemen- 
tary methods to identify bursts in the data, and it is thus in- 
teresting to compare the results from the two pipelines. We 
want to make the strongest statement that we can about the 
rate of gravitational wave bursts, consistent with sound statis- 
tical practice. We evaluated the efficiency of our two event 
trigger generator pipelines for each of our candidate signal 
waveforms, tuned to give approximately the same background 
rate using the playground data set. To combine the results 
from the two pipelines, we would choose to quote upper lim- 
its from the event trigger generator that has the better effi- 
ciency for the largest number of waveforms. With the tuning 
parameters used in this study, we find that the TFCLUSTERS 
pipeline has better efficiency than SLOPE for most of these 
waveforms (Table II), although SLOPE performs somewhat 
better for the 850 Hz and 1304 Hz sine-Gaussians. Thus inde- 
pendent of the decision to not use the SLOPE result because 
of the problems with background variability, the final upper 
limits that we would set are the ones obtained from the TF- 
CLUSTERS pipeline, shown in Fig. 15. The SLOPE results 
provide a valuable cross-check, and we intend to continue to 
employ and improve the SLOPE pipeline in future analyses 
(section VII B). 



C. Systematic uncertainties 

The following systematic uncertainties are known to sig- 
nificantly contribute to systematic errors associated with the 
results of our search. The estimation of our efficiency for de- 
tection of bursts with modeled waveforms (the Gaussians and 
sine-Gaussians that were considered in section V B) has a sta- 
tistical error associated with the finite number of simulations. 
The underlying noise floor was modeled using a sample of 
the S 1 data itself; there is a systematic uncertainty associated 
with the degree to which this sample was representative of the 
entire S 1 dataset. 

The analysis procedure for the efficiency has various po- 
tential sources of error. The detection efficiency versus h rss 
is fitted with a sigmoid curve, as discussed in section VB. 
The efficiency curves for each detector are combined to get 
an overall triple-coincidence efficiency, averaging over source 
direction and polarization, assuming that the coincident effi- 
ciency is the product of the individual efficiencies (as a func- 
tion of amplitudes at each detector). The effects of our (very 
limited) post-coincidence processing, including the choice of 
coincident time windows, clustering of multiple triggers from 



a single detector, and consistency of burst trigger parameters 
from the three detectors, have been studied. None of these 
studies reveal a significant source of error. It should be noted 
here that future results will employ far more detailed post- 
coincidence processing (see section VII B), including much 
tighter coincident time windows, and these issues will be care- 
fully re-evaluated. 

By far the largest source of systematic error associated with 
the efficiency determination is the uncertainty in the detector 
response to gravitational waves as obtained through the cali- 
bration procedure [27, 28]. We have combined the estimated 
uncertainty in the DC calibration with potential uncertainties 
in the frequency response, convolved with the frequency spec- 
tra of our modeled waveforms. We conservatively estimate an 
overall systematic uncertainty of 20% in the strain scale for 
our efficiency curves. Our curves of upper-limit versus grav- 
itational wave strain (Fig. 15) reflect this uncertainty by as- 
suming that the detectors' response is 20% less sensitive than 
the nominal value. 

We assign no systematic error associated with our choice of 
modeled waveforms since these are used to define the set of 
bursts which are targeted by this search. 

We assign no systematic errors associated with our choice 
of data used or our BLRMS cut or calibration cut (section 
II C). These procedures simply fix our observation time. 

The upper limit on the number of observed bursts depends 
on our estimate of the backgrounds, as discussed in section 
IV. The statistical errors associated with these estimations 
are marginalized over in the process of establishing the con- 
fidence bands; since these errors are small, this has a negligi- 
ble effect on the resulting upper limits. We have performed 
a variety of studies to search for and estimate sources of sys- 
tematic errors in the background estimate, including the time 
dependence of the background rate, various time-lag proce- 
dures, excluding anomalous data stretches, and other consis- 
tency checks. No sources of additional systematic error asso- 
ciated with the background rate are found to be significant. 



VII. SUMMARY AND DISCUSSION 

We have searched for gravitational wave bursts of astro- 
physical origin using data taken during the first Science Run 
of the three LIGO detectors. A total of 35.5 hours of triple co- 
incidence observation time satisfied our data quality require- 
ments. 

We employed two different algorithms for the identification 
of candidate burst events in the gravitational wave data chan- 
nel. With the algorithm for which we chose to present a final 
result (the TFCLUSTERS event trigger generator), we observe 
numbers of events that are reasonably consistent with expec- 
tations for random accidental coincidences of events originat- 
ing from noise fluctuations in the three detectors. We thus 
limit the excess event rate to be less than 1.6 per day , at 90% 
confidence level. 

We estimate our efficiency for the detection of linearly po- 
larized gravitational wave bursts incident on the detector ar- 
ray with a range of amplitudes and averaged over source di- 
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rection and wave polarization, by injecting simulated signals 
into the raw S 1 data streams and performing the search as it 
is done on the raw data itself. We focus on two simple, ad 
hoc waveform morphologies (section V A): Gaussians with 
a range of durations, and sine-Gaussians with a Q of 9, and 
a range of central frequencies. With these, we evaluate the 
(waveform-dependent) detection efficiencies as a function of 
strain per root hertz, h rss . We then combine the rate limit 
with the efficiencies to establish exclusion regions in the rate- 
versus-strength plane; representative examples are shown in 
Fig. 15. These constitute the results of the search reported in 
this paper. 

A. Comparison with other burst searches 



roughly h rss « 1 x 10 19 Hz 1 / 2 . To compare this with 
our sensitivity for 850 Hz sine-Gaussian bursts (h rss i/ 2 = 
7.3 x 10 -19 Hz -1 / 2 , Table II), we must first correct for our 
averaging over direction and polarization (section V C). This 
yields an amplitude at 50% efficiency for waves with optimal 
orientation of h rss w 2.6 x 10" 19 Hz~ 1/2 . Note that this 
detection threshold for LIGO is established by determining 
the loss of efficiency for fixed threshold, while for IGEC it 
is established by observing an increase in background events 
as the threshold is varied. Nonetheless, this measure of de- 
tection threshold permits a rough comparison of the search 
sensitivites [43], and we see that the IGEC search [4] has a 
somewhat greater sensitivity to 850 Hz sine-Gaussian bursts 
than the one presented here. 



Searches for gravitational wave bursts have been a main fo- 
cus of the observational program of the resonant-mass detec- 
tor community, and our work was able to benefit from their 
prior work. The most recent analysis by the IGEC consortium 
[4], which includes data from five detectors spread around the 
world, has presented its upper limits for bursts in the form of a 
rate-amplitude diagram, in much the same style as our Fig. 15. 

The IGEC result (Fig. 13 of [4]) bounds the rate for GW 
burst events with large amplitude to be less than ~ 4 x 10~ 3 
events per day. This rate bound is much stronger than the ones 
reported here due to the much longer observation time of the 
IGEC run. The rate bound gets rapidly worse for lower am- 
plitude bursts, due to increasing background and decreasing 
detection efficiency. 

To crudely compare the sensitivity of the two searches to 
low-amplitude bursts, we can define a "detection threshold" 
as the value of the burst amplitude at which the rate limit is 
two times worse than the limit for large amplitude bursts. For 
the analysis presented here, this is the 50% efficiency point 
h rS s 1/2 reported in Table II. 

One difference between the IGEC work and ours is that 
their instruments have relatively narrow frequency band- 
widths, and so are sensitive to a different measure of burst 
strength. Their detectors measure the Fourier magnitude 
\h(fb) \ of a signal waveform at the bars' resonant frequency 
fb ~ 900 Hz, incident at optimal source direction and polar- 
ization. 

However, if we consider a specific waveform with fre- 
quency content dominantly at or near at fb, such as our sine- 
Gaussians with central frequency /o = 850 Hz, the bar sen- 
sitivities and the interferometric detectors' sensitivities can be 
directly compared, over the relatively narrow frequency band 
where bar detectors are most sensitive. The conversion from 
| h(fb) | to h rss for bars for the sine-Gaussian family of signals 
(Eqn. 5.2) is calculated to be 



\kf b )\ 



ft (e^ 2 - r2 (f"-M 2 



(7.1) 



We focus on sine-Gaussians with central frequency /o = 
850 Hz, incident at optimal direction and polarization. Using 
Fig. 13 of [4] and Eqn. 7.1, the IGEC detection threshold is 



For all other waveforms shown in Fig. 15, and for other 
waveforms with significant spectral amplitude in a broad 
range of frequencies away from fb, the LIGO constraints are 
more stringent than the IGEC results, due to the broad band 
response of the interferometric detectors. For sine gaussians 
at 554 Hz, the ratio of peak spectral density of the pulse to 
the spectral density in the resonant-mass detector band is less 
than 10% for NIOBE and negligible in the other four IGEC 
detectors. For sine gaussians at 1304 Hz, ALLEGRO, AU- 
RIGA, EXPLORER and NAUTILUS receive spectral densi- 
ties in their bands that are only a few percent of peak spec- 
tral density for the pulse, with negligible spectral density in 
NIOBE's band. The resonant mass detectors also receive rel- 
atively small spectral density in their bands from gaussian 
waveforms compared to the LIGO detectors, unless r is less 
than 1 ms. This emphasises the importance of broad band sen- 
sitivity in searching for unmodeled gravitational-wave bursts. 
Ongoing work to broaden the response of resonant mass de- 
tectors should improve sensitivity to other waveforms in the 
future. 



The only previously published results on searches for burst 
events with broadband interferometric detectors that we are 
aware of are in Ref. [5] (but see also Ref. [6]). In Ref. [5], pro- 
totype detectors developed by the University of Glasgow and 
Max Planck Institute for Quantum Optics were operated for 
an effective coincident observing period of 62 hours in 1989. 
They searched for bursts with significant frequency content in 
the band from 800 to 1250 Hz. They considered the waveform 



h(t) 



; sin(27r f m t) j (2tt f m t) [44], which has constant 



Fourier magnitude from to f rn = 1250 Hz, and in-band 
(A/ = 1250 - 800 = 450 Hz) root-sum-square amplitude 
hrss = \/2A// (2f m )h pea k. They observe no events with 
h > h peak = 4.9 x lO" 16 , or h rss = 5.9 x 10~ 18 Hz~ 1/2 , 
averaging over wave polarizations and incident directions. 
Therefore, they set an upper limit on the rate of bursts with 
strain greater than this value, of 0.94/day. Their sensitivity 
can be compared with the strain sensitivities reported here at 
50% efficiency, for sine-Gaussians with central frequencies of 



850 Hz and 1304 Hz: h r 



i 1/2 



= 7.3 x 10 



-19 



Hz 



-1/2 



and 



l rss 1/2 



1.4 x 10~ 18 Hz~ 1/2 , respectively (Table II). 
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B. Directions for improved analysis in the future 

LIGO's second science run (S2) accumulated data for 8 
weeks in early 2003. At most frequencies, the noise in the 
three LIGO detectors was improved compared to the noise 
level of the SI data presented here by a factor of 10. Some 
improvements in the stability of the noise were also achieved. 
The in-lock duty cycles of the detectors were comparable to 
those obtained during SI, but tighter monitoring of the detec- 
tors' noise levels and calibration should lead to significantly 
less loss of data than was suffered in SI. Even without im- 
provements in our analysis methodology, we expect to ob- 
tain results from the S2 data that are an order of magnitude 
more sensitive in amplitude, and observation times that are 
increased by at least a factor of four over the results presented 
here. 

Based on lessons learned during the SI analysis, we are 
preparing numerous improvements and additions to our search 
methodology for the S2 data set. The pipeline presented here 
can be improved with more attention to optimizing and char- 
acterizing our event trigger generators. Obvious areas for 
improvement are better prefiltering, and better time resolu- 
tion for both SLOPE and TFCLUSTERS. We have also im- 
plemented an adaptive threshold for the SLOPE event trigger 
generator to make its event rate less sensitive to variations in 
detector noise. As the detector performance becomes more 
stable and closer to the design sensitivity, safe and effective 
vetoes based on auxiliary channels that monitor the environ- 
ment and interferometer sensing and control will be applied to 
reduce the number of spurious event triggers. New event trig- 
ger generators, using a variety of detection techniques (e.g., 
that proposed in [45]), will be tested. 

After initial (coarse) identification of coincident events with 
improved versions of the event trigger generators, the gravita- 
tional wave data channel time series can be reexamined to fur- 
ther reduce the background of accidental coincidences. Cross- 
correlation of the gravitational wave channels from multiple 
detectors can tighten our coincidence window so that it is lim- 



ited only by the light-travel time between detector sites, and 
test whether the event amplitudes and waveforms are consis- 
tent with the common origin of a gravitational plane wave. 

Future searches will include a more astrophysical style 
of interpretation, setting limits on populations of events in 
three-dimensional space. Efficiency simulations will include 
more realistic waveforms, such as black hole ringdowns or 
supernova waveforms [40, 41]. Higher sensitivity for mod- 
eled bursts can be obtained using matched filter techniques. 
Longer runs will give more opportunities to search for gravi- 
tational wave bursts coincident with gamma ray burst events, 
using the methodology described in [46]. 

Finally, and crucially, we are developing criteria by which 
we can establish confidence in the detection of gravitational 
wave bursts both statistically and as a single large amplitude 
burst event. For single burst event candidates, we will use 
information from all available detectors to reconstruct our best 
estimates of the gravitational wave direction, polarization, and 
waveform. 
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